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Abstract 

We consider (max,+)-algebra products of random matrices, which 
arise from performance evaluation of acyclic fork-join queueing net- 
works. A new algebraic technique to examine properties of the product 
and investigate its limiting behaviour is proposed based on an exten- 
sion of the standard matrix (max,+)-algebra by endowing it with the 
ordinary matrix addition as an external operation. As an application, 
we derive bounds on the (max,+)-algebra maximal Lyapunov exponent 
which can be considered as the cycle time of the networks. 
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1 Introduction 

We consider (max, +) -algebra products of random matrices arising from 
performance evaluation of acyclic fork-join queueing networks. The problem 
is to examine limiting behaviour of the product so as to evaluate its limiting 
matrix and the maximal Lyapunov exponent normally referred to as the 
system cycle time. 

In order to investigate the products, we develop a pure algebraic tech- 
nique similar to those involved in the conventional linear algebra. The tech- 
nique is based on an extension of the standard matrix (max, +) -algebra 
[H El El H] by endowing it with the ordinary matrix addition as an external 
operation. New properties of the extended algebra are then established in 
the form of inequalities, which may find their applications beyond of the 
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scope of the current topic. We conclude the paper with an example of ap- 
plication of the proposed technique to establish bounds on the cycle time 
and on its related limiting matrix in fork-join queueing networks. 

In fact, there exist similar results on evaluation of the Lyapunov expo- 
nent (see, e.g., [3] and references therein). However, they are essentially 
based on the description of system dynamics and related proofs made in 
terms of either Petri nets or stochastic events graphs. On the contrary, we 
exploit a different approach (see [5] for farther details) based on pure alge- 
braic techniques. It allows one to write and handle the dynamic equations 
directly without having recourse to an intermediate description in the Petri 
nets or in another tedious language. 

2 Motivating Example and Algebraic Model 

Consider a network of n nodes, with its topology described by an oriented 
acyclic graph. The nodes that have no predecessors are assumed to repre- 
sent an infinite external arrival stream of customers. Each node without 
successors is considered as an output node which releases customers from 
the network. 

Each node has a server and infinite buffer operating as a single-server 
queue under the first-come, first-served discipline. At the initial time, the 
servers and their buffers are assumed to be free of customers, except for 
the buffers in nodes with no predecessors, each assumed to have an infinite 
number of customers. 

The operation of each node can include join and fork operations which 
are performed respectively before and after service. The join operation is 
actually thought to cause each customer which comes into a node not to 
enter the queue but to wait until at least one customer from all preceding 
nodes arrives. Upon arrival, these customers are replaced by a new customer 
which joins the queue. 

The fork operation at a node is initiated every time the service of a 
customer is completed. It consists in replacing the customer by several new 
customers, each intended to go to one of the subsequent nodes. 

For the queue at node i, we denote the kth departure epochs by Xi(k) , 
and the kth service time by t^. We assume to be a given nonnegative 
random variable (r.v.) for all i = 1, . . . ,n, and = 1,2,... 

We are interested in evaluating the limit 

1 

7= lim -maxxj(fc), 
k— >oo k i 

which is normally referred to as the cycle time of the network. 

In order to represent the network dynamics in a form suitable for fur- 
ther analysis, we exploit the idempotent (max, +) -algebra based approach 
developed in [5]. 



2 



The (max, +) -algebra [H El [3] presents a triple (i? £ ,ffi, <g>) with R e = 
RU {e}, e = — oo, and operations and <g> defined for all i,y£jJ E as 

x © y = max(x, y) , x ® y = x + y. 

The (max, +) -algebra of matrices is introduced in the ordinary way. 
The square matrix £ with all its elements equal e presents the null matrix, 
whereas the matrix E = diag(0, . . . , 0) with e as its off-diagonal elements 
is the identity. 

Let us denote the vector of the feth customer departures from the net- 
work nodes by x(k) = (x\(k), . . . ,x n (k)) T , and introduce the matrix Tk = 
diag(Tifc, . . . ,r nk ) with all its off-diagonal elements equal e. 

As it has been shown in [5], the dynamics of acyclic fork-join networks 
can be described by the stochastic difference equation 

v 

x(k) = A(k)®x(k-l), A(k) = Q){T k fg> G T ) j <g> T k , (1) 

3=0 

where G is a matrix with the elements 

J 0, if there exists arc in the network graph, 
^ \ e, otherwise, 

and p is the length of the longest path in the graph. 

The matrix G is normally referred to as the support matrix of the net- 
work. Note that since the network graph is acyclic, we have G q = £ for all 
q>p. 

Consider the service cycle time 7. Now we can represent it as 

7 = lim d|as(&)||, 

k— >oo K 

where ||s(/c)|| =maxjXj(/c). 

In order to get information about the growth rate of x(k), we will ex- 
amine the limiting behaviour of the matrix 

k p 

A k = A(k) ® • • • ® A(l) = (g) 0(T fc ® G T y ® T k . 

i=i 3=0 

3 Distributivity Properties and Matrix Products 

Let Aij be (n x n)-matrices for all i = 1, . . . , k and j = 1, . . . ,m. Distribu- 
tivity of the operation ® over © immediately gives the equality 

k m 

(8>0^= © A ln ®-..®A kjk , (2) 

i=l j=l l<ii,-..j'fe<m 
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which leads, in particular, to the inequality 

km m k 

(g)©^>©(gM*r ( 3 ) 

i=l j=l j=l i=l 

We consider the ordinary matrix addition + as an external operation, 
and assume (g> and © to take precedence over + . In a similar way as above, 
we have 

k m 

E0^i = © (A ln + ... + A kjk ), (4) 

i=l j=l l<ji,---,jk<m 
km m k 

/s&Ajj > G^/.Ajj- ( 5 ) 

i=l j = l j = l i=l 

Let Gi and G 2 be support matrices. For any matrices A and B, we 
have 

G 1 ®(A + B) ®G 2 <G 1 ®A®G 2 + G 1 ®B®G 2 . (6) 

Assume Di and D 2 to be diagonal matrices with all off-diagonal ele- 
ments equal e. Then for any matrices A and B , it holds 

D X ®(A+B)®D 2 = Di®A®D 2 + B = D X ®A+B®D 2 = A + Di®B®D 2 . 

(7) 

Now we examine products of alternating diagonal and support matrices 
denoted respectively by D and G, which take the form 

D <g> (G ® D) ® • • • ® (G ® -D) = £> ® (G ® D) m . 

V v ' 

m times 

In order to simplify further formulas, we introduce the following nota- 
tions 

= D ® (G ® £>) J ', = G* <g> -D <g> G J . 

First assume the diagonal matrices to have both positive and negative 
entries on the diagonal. The next lemma can be proved using ([6]) and in- 
duction on m. 

Lemma 1. It holds that 

m 

<S>m(D)<Y,*7~ 3 (D). 

3=0 

Furthermore, assuming , i = 1, . . . , A;, to be diagonal matrices, one 
can obtain the next result based on Lemma [T] and inequality dU]). 
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Lemma 2. Let mi, . . . , writ be integers, and m = mi + • • • + m^ . Then it 
holds 

k k Mi 

(g)$ mi (A)<E E *T s &i) 

i=l i=l j=Mi-i 

with Mq = , Mi = mi + • • • + rrij, i = 1, . . . ,k . 

Let the matrices Di , . . . , A have only nonnegative elements on the di- 
agonal. With ([7]) and ([6]), one can prove the next lemma. 

Lemma 3. Suppose that mi + • • • + m r = m r +i + • • • + m^ = m with 
m — m r < m r+ i for some r . Then for any integer s such that m — m r < 
s < m r+ i , it holds 

r k k 

(g)$ mi (A) + (g) $ mi (A)>(g)^(A) 

i=l i=r+l i=l 

with si + ■ ■ ■ + Sk = m, and 

mi, ifl<i<r, 

s — m + m r , ifi = r, 

m r+ i — s, if i = r + 1, 

mi, if r + 1 < i < k. 

4 Subadditivity Property and Algebraic Bounds 

Consider the family {Af k \l, k = 0, 1, . . . ; I < k} of matrices 

Aj k = A T (l + 1) ® • • • ® A T (k), A T ok = Al 

The next lemma states that the family {Af k } possesses subadditivity 
property. 

Lemma 4. For all I < r < k, it holds 

Afk < Af r + Af k . 
Proof. By applying ([2]) and @ , and then Lemma O we have 

r p k p 

Af r + A? k = (g) 0^(73+ (g) ©*i(7i) 

i=l+l j=0 i=r+l j=0 

> © <g> 

m=0 s;+iH \-s k =m i=l+l 

Finally, since G m = £ for all m > p, we get 

k k p 

Al + A T rk > (g) $ Si (%) = (g) $ s {Ti) = Al. □ 

0<s l+l ,...,s k <pi=l+l i=l+l j=0 



5 



The next lemma offers bounds on . 
Lemma 5. It holds that 

[p/kj k k p k 

©C8)^(^)<^<||07i||®0G r +^ Km), 

r=0 i=l i=l r=l i=l 0<r+s<p 

where \r\ denotes the greatest integer equal to or less than r . 

Proof. The lower bound is an immediate consequence from ([3]), and the 
condition that G m = £ if m = kr > p. 

In order to derive the upper bound, we first apply ([2D to write 

k 

A k = ®*m*(70- 

0<mi,...,m,fe<p i=l 

Application of Lemma [2] gives 

k k M- 

®wti)<i; e *r s w 

i=l i=l j=Mi-i 

k M: k 

= E E *r j ra + £*M-f-m 

i=l j=Mi_i+l i=l 

With ([5]), we further obtain 

fe M, fe 

^< EE *r i (7i)+ © w 

0<mi,...,m,t<p i=l j=Mi— i+l 0<mi,...,m^<p i=l 

It remains to replace the first sum with its obvious upper bound, and 
then apply and ([5]) to the second sum so as to get the desired result. □ 

5 Evaluation of Bounds on the Cycle Time 

The next statement follows from the classical result in [BJ, combined with 
Lemma 01 

Theorem 6. If Tn,Ti2, ■ ■ ■ , are i.i.d. r.v.'s with E[rji] < oo for each i = 
1, . . . ,n, then there exists a fixed matrix A such that with probability 1, 

lim Al/k = A T , and lim E[A%]/k = A T . 

k— >oo k—^oo 

Furthermore, application of Lemma [5] together with asymptotic results 
in [8] gives us the next theorem. 
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Theorem 7. If in addition to the conditions of Theorem® ~D[th] < oo for 
each i = 1, . . . , n, then it holds 



E[71] < A T < E 



G r ®Ti®G s 

. 0<r+s<p 

As a consequence, we have the next lemma. 



(8) 



Lemma 8. Under the conditions of Theorem^ for any finite vector x(0) , 
it holds 



|E[Ti]|| <7< 



E 



G T ®Ti®G s 

0<r+s<p 
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